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I. Introduction 

In a recent paper [1] (hereafter referred to as I) we have obtained the RPA equations for 
0^ field theory by hnearizing the time-dependent variational equations. The method was 
implemented for the case of symmetric vacuum, (0) = 0, which allows us to investigate two- 
meson physics. We have shown that it is a simple nonperturbative method to study scattering 
processes. Using this framework the problem of stability of vacuum can be explored from the 
RPA modes. In continuation of I we will consider here the stability of the theory for other 
critical points of the Gaussian parameter space [2]. In particular, we discuss the solutions 
and renormalization conditions for the asymmetric vacuum. In this case, the excitations of 
the vacuum are identified with the one and two-particle wavefunctions and the system of 
equations can also be solved analytically. As result, a stable zero mode is found for certain 
range of renormalized coupling constant and as well as a complete form of the scattering 
amplitude. 

For completeness, let us first repeate here the key equations of I. The bare parameter 
hamiltonian for the 0^ theory is [We use the notation: /x = / d^x] 

H = lil^'i^) + livmr + |0^(x) + ^0^(x)) . (1.1) 

As an approximation, we take a Gaussian trial wave functional 

il;{cl>,t) ^ Nexpl - £ [^(x) - 0o(x, t)] f ^' - zS(x, y, t)] [^(y) - 0o(y, t) 



+ i£7ro(x,t)[0(x) -0o(x,i)]l, (1.2) 



where N is a, normalization and 0o(x, t), 7ro(x, t), ^(x, y, t) and E(x, y, t) are our variational 



parameters. These quantities are related to the following mean-values: 

(V^,t|0(x) I V^,t) = 0o(x,t), 
I 7r(x) \^,t) ^ 7ro(x,i) , 
t I 0(x)0(y) I V, ^) = ^(x, y, t) + M^, t)My, t) , 
{^,t I I = / E(x,y,t)G(x,y,t) + [ 7ro(x, t)0o(x, t) 
From (1.1) and (1.2) one can compute the effective hamiltonian, 



(1.3) 

(1.4) 
(1.5) 
(1.6) 



n = {i^,t\H\^,t), 



(1.7) 



which is the energy of the system, 

n = £[^7r2(x,i) + ^(V0o(x,i))' + |0^(x,i) + A04(^^^y 

+ i / ^-^(x, X, t) + 2 / E(z, X, t)G'(x, y, t)S(y, z, t) 

+ ^ £[-V^ + a + ^(/.2(x, t)]G'(x, y, t) + ^ £ ^(x, x, t)G(x, x, t) . 

The variational equations of motion read 

0o(x,i) = 7ro(x,i), 

7ro(x,t) = VVo(x,t) - a0o(x,t) - -0o(x,t) - -0o(x,t)G(x,x,t) , 
G'(x,y,t) = 2 /fG(x,z,t)S(z,y,t) + S(x,z,t)G'(z,y,t) 
E(x,y,f) = -2^E(y,z,i)E(z,x,i) + ^G'-2(x,y,i) 



(1.8) 



(1.9) 
(1.10) 
(1.11) 



+ 2 



- a - -</.^(x, i) - -G(x, X, t)] 5(x - y) . 



(1.12) 



These are nonlinear time-dependent field equations. Therefore, a closed solution is not 
easily constructed. Here, we will consider the equations in the equilibrium situation and the 
small oscillation regime. In these cases, an explicit solution can be obtained allowing us to 
examine diverse properties of the theory. The structure of this paper is as follows. Section 
II discusses the time independent Gaussian equations. Particular attention is paid to the 
question of the renormalization and we study solutions for the resulting gap equation. We 
also discuss the so called Gaussian effective potential for the case when the bare couphng 
constant 6 — > 0"*". Section III discusses stability of the stationary points found in the section 
II by investigating the properties of variational space. In Section IV we shall derive the 
RPA equation by considering near equilibrium dynamics about the critical points. Since the 
excitations of vacuum are quantum particles. Section V will solve the RPA equation with 
boundary conditions of scattering processes. In Sec VI we will use the properties of separable 
potential to get solutions for the T matrix. The spectrum of RPA modes is discussed using 
this result and conditions of stability will be analysed within this context. 



This section will discuss the variational equations, Eqs.(1.8)-(1.12), in the equilibrium 
situation. We investigate the solutions of these equations and study the renormalization 
conditions. We will consider 4>q uniform, 4>o{x) = ip, and parametrize the kernel G as 



independent of x. With these assumptions the expectation value of H, (1.8), with kinetic 
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II. Time-Independent Gaussian Variational Equations 




Therefore, the parameters ip and m are taken to be 



term set to zero, reads as 

U{(p,m) 
u{(fi, m) 



m] 



1 ?T7^ 1 r h i2 

-G-^(m)-— G(m) + -[a + -G(m)^ 



1 



b 
2* 



(2.1) 



In this equation G{m) and G ^(m) are defined as (x | G | x) and (x | G ^ | x) respectively. 



By virtue of parametrization of G these matrix elements are given by 



G{m) = (x 



1 



x) 



1 



G ^(m) = (x I 2Jp'^ + m? \ x 



1 



kl-mHog' ^"^P 



87r2 



\fem' 
2K 



^ ^ Vem 4 

where Ap is a cutoff in the integrals. Notice that the energy density m is a function of two 
variables </? and m. The next step is to obtain the critical points of this functions and find 
the renormalization conditions. 

II— a. Gap Equation and Renormalization 



(2.2) 
(2.3) 



Minimization of u with respect to </? and m yields 



c)Vb r h b 

— — = mG'{m) a H — H — G{m) — w? = , 

dm '-2 2 J 



(2.4) 
(2.5) 



where G'{m) = dG{m)/dm'. Notice that one may have several sets of solutions for the 



system of equations (2.4-2.5). 



By analysing the stability matrix (see section III) one can show that the solution m = 
of (2.5) is always a saddle point of u{ip, m) because the determinant is negative in this case. 
Thus, the possible stable points of (2.1) are determined through the gap equation, defined 
by 

m' = a+^^' + ^Gim), (2.6) 
and (2.4). From its solutions, (f = and 

a + ^^' + ^G{m) = 0, (2.7) 

one can eliminate ip in favor of m and arrive at the following key equations to be solved, 

= a(a + ^G(m)) . (2.8) 



In this expression, we have parametrized the two solutions of (2.4) as follows: a = 1 when 
(f = and a = —2 for b(p'^ = —6[a + b/2G{m)]. On the other hand, by combining ( p.6|) and 
( p.7| ) it is easy to show that 

bip^ = (2.9) 

which we define as the broken phase. Notice that this phase might degenerate with the 
massless symmetric phase if m = is a solution of (T^). We will come back to this point 
later. 

Let us now examinate the renormalization conditions for the gap equation (p.8|) . First, 
we use (p. 21) to rewrite (WM as 



bAi bm^ , 2Ap 

— = « + - TTT^ log -f^ ■ 2.10 
a IStc^ Idtt^ ^/em 



Next we introduce an arbitrary mass scale /i by writing log(Ap/m) = log(Ap/yu) +log(/i/m) 
the equation (|2.1(]| ) becomes 



= 1 \^^^%A + J f 2A ^ • (2.11) 

This equation renders finite result if we choose the following self-consistency condition 

h\2 



= 1 , '""^ OA ^ e = 0,±l, (2.12) 



a+ P 
IGtt"^ 


1 




2Ap ' 




^ 1 « 2 




a 


lOTT^ 


v/e/i 


2 

a6 







12 1 2Ap , , 

log^. (2.13) 



ao ■ ■ avr^ " e/i 



Thus, the renormalized version of Eq. (p. 81) reads 



— 2 

= e/i^ + ^^^m^ log ^ , (2.14) 

where = (7^/167r^ and e indicates different renormalization conditions. The above result 
can be interpreted as follows. We start with two bare parameters a and b. With the help of 
the transformations rules (|2.12|) - (p.l3|) we arrive at a new set of paramteters fi and g. Being 
fi a renormalized mass, which also defines a mass scale, and g a dimensionless renormalized 
coupling constant, in a such way that the resulting gap equation involves finite quantities 
only. Notice also from (|2.12|) - (|2.14|) that two different theories of (<^^)3+i are involved. One 



is the version with 6 oc 1/L^ (a = 1) and the other is 6 oc — 1/L^ (a = —2) which belong 
to two distinct field theories [3-4] . In the limit of infinite cutoff the renormalized coupling 
constant g^ can assume any value. 

II— b. Massive and Massless Solution 
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In order to discuss the possible solutions of (2.8) and investigate the roles played by e we 
compute u{ip{m),'m), being iflm) the solutions of (2.4), 

M(m) = u{ip{m),m) 

G'\m) - —G{m) - + -G{m)] (2.15) 



4 ' ' 2 ' ' 26 

1 / 4i 4 , 

= ^ ^ log — : — m + jL , , , ^ , 

1287r2V ^ ^4 y ^ Q^tt^ 

where we have used (2.3) and the counterterms introduced in (2.12)-(2.13)[1,2]. The different 
behavior of u{m) with < and e = 1 are shown in the fig. la of I. In this case u{m = fi) 
and u{m = 0) are local or true minima depeding on the values of However m = shown 
there actually corresponds to the massless solution of (|2.5|) , which is a saddle point in the 
{(p,m) space. In addition, u{m) calculated with (7^ > reduce to those in fig. la if one uses 
the scale rh defined as minimum of u{m?) (see appendix B of I) . A similar consideration 
can be made also for the case e = — 1 as follows. 

Recalling ( |2.13| ) we can write a relation of renormalized coupling constants at mass scales 
/i and fh readly as 

1 1 , . X 

— = --log — , (2.17) 



where m is defined as the minimum of u{m). Combining now (p.l4|) with (p.l7|) one can get 
the following simple relation 



(l.i-)^V(l.l). (2.18) 
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Using this result we can rewrite (|2.16|) , at scale of m, as 



u[m) = — log —---— + 1 — - if. 2.19 



Disregarding some unimportant additive constants this equation is equal to (|2.16| ) in unit of 
m. Therefore, this renormalization scheme does not add any new physics. 

The case when e = corresponds to the renormalization prescription a + bA'^/lGir^ = 
[2,5]. Thus, the gap equation (2.8) becomes 

^'(- + T7^1og4^)=0 (2.20) 



which allows a massless solution. It is important to notice that there is only one free param- 
eter involved in this case. Thus, we may fix it using the mass scale rh defined by the second 
solution of (|2.20|) , namely 

ab 877^ ^ye^m 

Comparing this to ( p. 131) and using (|2.17|) one has (jm = —1- With these ingredients we find 



the following the renormalized version of u for this case: 



4 4 

TTl Tf) 



On the other hand, one can obtain the above result rewriting directly ( |2.16| ) in unit of 



defined by ( p.20|) . In other words, the case of e = corresponds to the result given by 
( p.l9| ) with = —1- Therefore, we can conclude after this analysis that all possible physics 
situations in u{(p{m),m) are contained in the case of Qm < and e = 1 when rh is utilized 
as the mass scale. It has minima at m = m, and massless solution m = 0, for the particular 



case of = —1. Thus, the time-independent variational equations ( p.4|) -( p3D allow simple 



analytic solution and energy density u{ip, m) has equilibrium points at {ip = 0, m) and 
{bip"^ = 3m^, m). 

II— c. Gaussian Effective Potential 



In the previous subsections we have studied the energy density u{(p, m) when it is min- 
imum with respect to ip, i.e., u{ip{'m),'m), being ip{m) given by (2.4). Now we want to 
calculate u{ip,m) when it is minimum with respect to m. In particular, when the resulting 
expression is written as a function of cp, it is known as Gaussian Effective Potential (GEP) 



[3-4]. However, an explicit expression of m((p) from ( |2.6| ) is not straighforward. Instead, we 
write ip in terms of m, 

hn^ „ h 

(2.23) 



V 2 b 

= m —a G m . 

2 2 ^ ^ 



Substituting this into (pTT| ) yields 



-G{m) 



Dim) 



(2.24) 



where 



Dim) 



m^ 2m^ 



66 ^ 36 



a + -G(m) 



+ 



36 



a + -G(m) 



(2.25) 



Notice that u{m) — D{m) is exactly equal to ( p.l5[ ) for the of case a = —2. Using now 
the renormalization procedure given by ( p.l2|) -( p.l3D it is straightforward to show that D 
0(^-^^. Therefore, the energy density at the curves defined by ( ^l6|) as well as by ( p.7|) in the 
((/?, m) plane has the same result. The coincidence is not casual and let us look at (2.1) more 
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carefully. Recalling the relation (|2.9|) we note that since is a finite physical quantity, 
therefore it suggests that (p requires a scaling factor. Thus, we define a finite mean-field 



value as 



(2.26) 



Using now ( p.l2| )-( p.l3| ) one can rewrite ( |2.6| ) and ( p.7| ) respectively as 

1 



$2 = ^ 



$2 — 0™2 



rn? — ij? 



3m^ + 



167r2 



16vr2 V^i^ 



m m 

, mi 
^log- 



(2.27) 
(2.28) 



This result shows that in the limit of infinite cutoff, the two curves converge to bip"^ = Sm^. 



Using this relation one can write GEP from ( p.l6| ) immediately as 



$4 $2 

5767r2 V °S ^ 



1 N $2 .$2 

^^^^ 1927r2VY 



2/i^ . 



(2.29) 



In particular, (2.29) recovers the GEP obtained by the authors of Ref. [3] after some changes 
of variables (cf. their Eq.(17)). 

III. Stability Analysis 



In this section we shall analyse the stability conditions for the solutions obtained in the 
previous discussion. Since the energy density u{ip,m), (|2.1| ), is a function of two variables, 
we can define its stability matrix as 



A 



difdif dipdm 

\ dmdip dmdm 
a + ^G{m) + |(/?2 
bipmG' {m) 



bipmG'im) 

[G'{m) + 2m2G"(m)]r + 2m^G'{ni)[\G'{m) 
11 



1] 



,(3.1) 



where 

V = a+^-if^+^-G{m)-m\ (3.2) 

Now the stabihty analysis reduce to a discussion about the signs of the eigenvalues of A at 

u{ip, rh), being and rh the solutions of (^.4|)-(P3|). 

Section II has showed that when F = 0, the system might have two solutions at y9 = 

and = 3m^. For the former case A is diagonal, 

/ m2 \ 

^ " I, 2m''G'{m)[^G'{m)-l] )■ ^^'^^ 

Hence ip and m are the eigenvectors of A. In addition, is the oscillation mode in the (p 
direction. Recalling now the renormalization condition ( |2.12| )- (|2.13|) with a = 1 we get the 



following expression for the eigenvalue in the m direction, 



The system is stable if 

1 1 



<^log — . (3.5) 



This is the result we have obtained in I from the RPA analysis. 

For the broken phase the stability matrix is no longer diagonal, 

/ b(prhG'{rh) \ 

^ ^ [biffhG'im) 2m^G'{m)[^G'{m) -1] ) ' ^^'^^ 

Its determinant is equal to 

detA = -2m'^G'(m)(bG'(m) + l) = 2mH-— + ^log^), (3.7) 
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where we have used the relation h>^'^ = 3m and the renormahzation condition ( p.l2|) -( p.l3|) 



with a = —2 in the second equahty. The criterion of minimum requires the eigenvalues to 
be positive. A straighforward calculation yields 



A+ = 



2m' 

A 



L,{3 + 0{M}, (3.8) 

;r + T;^;;i'°8TT) + o(7j)- (3-9) 



The eigenvalue A"^ is always positive while A^ > if the Eq.(3.5) is satisfied. The result 
shows that u{(p,m) is quite singular at this critical point, because A~ is a zero mode for 
the limit of infinite cutoff and the curvature correspondent to the eigenvalue A"*" is infinitely 
sharp. 

IV. RPA Equations 

Section II discussed in detail the possible minima of 7i as a function of the variational pa- 
rameters. In this and next sections we will investigate near equilibrium dynamics around the 
stationary points. In I we have obtained the RPA equations by linearizing the time-dependent 
variational equations. Here we will procede differently: we first expand the hamiltonian, Ti, 
around the stationary points and then the RPA equations are obtained from Hamilton's 
equation with the new hamiltonian [6]. Of course these two approachs are equivalent, but 
the method discussed here has advantage of eliminating the step of introducing a new aux- 
iliary variable 6v (see Section IV of I). 

Let us first consider fluctuations around <^ and rh, which are generic solutions obtained 
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in the previous section for the uniform system, 



5G{^,y,t) = ^(x, y, - G(x, y) . 



(4.1) 
(4.2) 



They define the one- and two-particle wavefunctions respectively. These quantities and their 
canonical conjugate momenta, ttq and E, are assumed to be small in our approximation. 
Next, we expand the hamiltonian, Eq.(1.8), up to second order in S4>, ttq, SG and E, 



n 



RPA 



+ 
+ 



'\ [a + ^G(x, x)]<^^ + (x I \a~" + \fG I x) + 1 [a + ^G(x, x)]^ + A^^^ 
a+^(p^ + ^<5(x, x)](^ j (50(x, X, t) 



/,(_ 

Ly (2 

+ \j 7r2(x, t) + 2 / E(z, X, t)G(x, y)E(y, z, t) 

I J:x. •'XjYjZ 



- + a + \(p^ + ^G(x, x)] 5(x - y) - ^^-^(x, y)] 5G(x, y, i) 



+ 



h ^ 



hip 



-V^-l-a-t--<^' + -G'(x,x)j50^(x,t) + :^ / 50(x,t)5G(x,x,t) 



+ ^ / 5G'(x,x,t)5G(x,x,t) + - [{x\G \gG \gG ^ | x) . 
8 >/x 8 7x 



(4.3) 



The first term of this expression is the time-independent part of the hamiltonian. The second 
and third terms vanish because of the equilibrium conditions, (2.4) and (2.5). The Eq. (2.6) 
can also be used to simplify the sixth term. The matrix element of the last term can be 
written explicitly in momentum space as 



X I G ^5GG \gG ^ I x) 



= / UJp^UJp^{uJp^+UJp2)^G{pi,P2,t)SG{-pi,-p2,t) 



(4.4) 
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where cUp is a matrix element of the operator in momentum space, 



(pi I a) I -P2) = (pi I ^/p^iPm^ I -P2) = 5{pi + P2)a;p. 

It is useful to introduce the relative and total momentum coordinates given respectively by 

p = (pi-p2)/2, (4.5) 
P = (P1 + P2). (4.6) 

Thus, from these remarks the Eq.(4.3) can be rewritten, in terms of p, P coordinates, as 

Hrpa = ?iE5T+^ / 7ro(P,i)7ro(-P,t)+ / E(p,P,t)(-^ + -^)E(-p,-P,i) 

S<P{P, t)5G{p, -P, t) + 5(f){-P, t)5G{p, P, t) 



i 


1 


p2 




2 




hip 






T 


J 


p,P 


h 


[ 




8 J 


P,P',P 



+ ]-[ a;+a;_(cu++a;_)5G(p,P,t)5G(-p,-P,t), 
2 Jp,p 



(4.7) 

where u}± = y^(p ± P)^ + m^. It is convenient to make the following changes of variables for 
the further use. 



E(p,P,i)/i(p,P,i) = n(p,p,t), 
5G{p,P,t)h-\p,P,t) = p{p,P,t), 



(4.8) 
(4.9) 



being 



h{p,P) 



(4.10) 
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Thus, 



Hrpa = Hest + \ f M'P,t)no{-P,t) + l [ n(p,P,i)n(-p,-P,i) 
2 Jp Z Jp,F 



+ 
+ 



[- 


P 


Jp 2 








4 J 


P,P 



h{p, P) [50(P, t)5p{p, -P, t) + 5cj>{-P, t)5p{p, P, t) 

+ U (h{p, P)Sp{p, P, t)) (hip', P)dp{p', -P, t 

+ i / (a;+ + uJ- fSp{p, P, t)5p(-p, -P, t) . 
2 yp,p 



(4.11) 



From this one can get the hnearized equations of motion by following directly Hamiltion's 
equations [1]: 



7ro(P,t) 
Sp{p,P,t) 

ri(p,p,i) 



- P^ + rh^ 

n(-p, -p,i), 



6(P{-p, - Y Hp, P)6p{p, -p, t) , 



(4.12) 
(4.13) 
(4.14) 



+ ou^)'5pi-p, -P, t) - ^-h{p, P) h(p', P)5p{p', -P, t) 
f h{p,P)54>{-P,t). (4.15) 



b(p 



Eliminating the canonical momenta ttq and U we arrive finally at 



hip 



5(/.(P, t) + ul.6<P{P, t) + ^ J^^ h{p', P)Sp{p', P,t) = 0, 



(4.16) 



Sp+ {uj++uj.f6p{p, P, t) + -h{p, P) h{p', P)6p{p', P, t) + -^h{p, P)6<P{P, t) = 0. (4.17) 

These are linear oscillator equations as usual in RPA treatment. The solutions for this 
problem involves determining the modes of oscillation and their eigenfrequences. 
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We have seen in Sec. II that (p may have solutions at (y9 = and <y9^ = "iw? jh. Note that 
in the symmetric phase the two equations decouple. The two-particle equation reduces to 
Eq.(4.20) of I while the one-particle equation becomes simple oscillator equations, one for 
each P with frequence ti;p, 

50(P, t) + cu|50(P, t) = . (4.18) 

For the solution Lp^ = Srh'^/b, however, the coupling between the one-particle and two-particle 
wavefunction is nontrivial, we shall discuss the solution of this problem in the next section. 

V. Small Oscillation Equation and Scattering Problem 

Let us first remove the time dependence of ([4.16| )-( P:.17| ) by writting 



S(j){P,t) = 50(°)(P)cos[fi(P)t], (5.1) 
6p{p,P,t) = 5p(°)(p,P)cos[l](P)t], (5.2) 

where 5(f)^^^ and Sp^^^ are the amplitude of oscillation for the one- and two-particle modes 
respectively and ^^(P) are the eigenfrequences (see section IV-b of I for interpretation of 
the wavefunctions) . Substituting these into ( |4.16|) - ([4.17|) we have the following eigenvalue 
problem: 

^'M'^ (P) + Y l Mp, p)¥'^ (p', p) = ^'H^"^ (P) , (5.3) 



[UJ- 



+ a;_)V°Hp,P) + h{p,P)l h{p',P)Sp^'\p',P) 



+ ^/i(p,P)50W(P)=f]V°nP,P)- (5.4) 
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Since the amplitudes of oscillations are wavefunctions of quantum particles it is interesting 
to treat this system as a coupled channel scattering problem with appropriate boundary 
conditions. Henceforth, we shall use a and /5 to denote the one and two-particle channel 
respectively. In the following discussion we will consider the cases separately. 

a. a — *■ a 

In this process one has a incident wave of 50^°-'; it couples to 5p'^°^ through the term 
hip J hSp^'^^ and reemits 50'^°^ at the exit channel. Thus, we can formally solve ( |5.4] ) as 
follows: 

V°nP,P) = o2 ( ]_ ^2^- [iMp.P) / h{p',P)Sp{p',P,t) + ^Jh{p,P)6<p{P,t) . 

ii^ — (ci;_|_ + C(;_j^ + '-4 Jp / J 

(5.5) 

This expression includes the boundary condition that there is no indent wave of Sp^'^\ The 
term Jph{p,P)5p^^\p,P) couples 5p^^^ of different relative momenta. We can solve for this 
term by multiplying (|5.5| ) by h{p, P) and integrates it with respect to p, 

/^Mp',P)<5p(°)(p,P) = ^^^^50(°)(P), (5.6) 



where 



The details of the computation of this integral can be found in the appendix B of I. For the 
present analysis it is sufficient to say 

/p(fi) = -L^ + Fp(fi), (5.8) 



where is the logarithmic divergent term defined in (|2.12|) and Fp{Q) the finite part of 



/p(r2) [See (6.3)-(6.6) for /p(i7)]. Substituing now ( |5.6|) into (|5.3| ) we have 



2/b-ip{n) 



50(°)(P) = 0. 



(5.9) 



For the asymmetric solution, bip"^ = 3m^, and ([5 .91) becomes 



1]2 _ p2 _ ^2 



(1 



2/b-Ip{n) 



50(0) (P) = 0, 



(5.10) 



Using the renormalization condition ( f^.l3D and ( |5.8| ) one arrives at 



P) = 0, 



(5.11) 



Notice that the effect of the coupling 5(j)^'^'> to the two-particle modes switch the mass of the 
particle. In the limit of infinity cutoff its effective mass goes to zero. 



b. /3 — > 13 

It is also an elastic scattering process, where the entrance channel as well as the exit 
channel exhibit two-particle wave ^p'-''-*. Analogously to the previous discussion we first solve 



(5.3) as 



(50(°)(P) 



b^/2 



Mp',P)5p(p',P). 



(5.12) 



fi^ — cUp + ie Jp' 

Using this and (5.4) we get the following integral equation for 5p^^^: 

jp'°'(k,p.P;n) ^7g<(k-p)+j( „, ^"^tl' ■ )''(P.P) / A(p',P)M°'(k,p',P;Si) 



(5.13) 
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In this expression k is the relative momentum of the two incident meson; 7^ indicates a 
phase factor to keep ^p'-*'-* real and the subindex /5 denotes the source of incident wave. From 
( ^.13| ) one finds immediately 

h{p, P)5pW (p, P) = 7^^^J— /,(p, P) , 



where 



Using this result in ( ^.13| ) we find 



5pW (k, p, P; n) = 7^5(k - p) + ] ^^j hik, P)S(P, n)h{p, P) (5.15) 



with 

Of course, this result reduce to Eq.(5.5) of I in the case of = 0. For the broken solution, 
we can rewrite ( [5.16| ) in terms of renormalized parameters as 



5(P,fi) = ^ . (5.17) 



3L^(r]2 - P2) - —iSf - P' - m^) - F(P, - P2 + 2m2] 



Notice that S L^^ for large Ap, which means that the scattering waves is infinitesimaly 
small. Therefore, effectively there is no interaction between the particles coming from this 
broken vacuum. 

c. a — > j3 
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In this case we want to obtain 5p^^^ with a source of one particle, which is an inelastic 



processes. To do so, (|5.3|) is solved as 



i(°)(P,K,,fi) =7,5(K-P) + 



/ h{p',P)6p^'\p',P,n). (5.18) 



Here we have included an incident wave of 50^°^ with moment K. Using this solution in (|5.4| ) 
one gets an integral equation for 5p^^^ with a source in a, 



7«<5(K-P)^Mp',P) 



+-(l + ^ 



Mp,P) / /i(p',P)5pW(p',P)(5.19) 



This equation can be solved as usual given 



5pW(K,p,P;fi) = ^J{K-P)^h{p',P)[l + E{P,n) 



7„5(K-P)^Mp',P) 



(5.20) 



In the second line we have used ( ^.ITj) . Recalling the discussion in the section II, we have 



learned that in the broken solution, the mean-field (f requires a scaling factor. As conse- 



quence, we have essencially b(^ oc l/JL^, which vanishes for large cutoff. Therefore, in the 



continuum limit 6p cannot be observed in this reaction. 



d. P — > a 



For the last case, the Equation (|5^ ) is solved with the source in (3, 



Sp^'\k,p,p-n) = 7^5(k-p) + 



1 



^hip, P) J^^ h{p\ P)5p(°)(k, p', P; + ^/i(p, P) 



(5.21) 
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From this, one gets immediately 

|^/i(p,P)5pW(p,P) = 1 _ 5jp(^)/2 h^^'^) + b^H^)S<PP^'\p)] ■ (5-22) 
Combining this with { f).3\j results 

^^"Hp) = iMHk) ^ J^^^,^ . (5.23) 



2 _ 2 _ b^^ipjn) 



Above equation differs ( ^.91) from the boundary condition for this scattering process. Here 
we have incoming waves of 6p characterized by 7^ in its rhs. It is proportional to bi^, or 
50*^°^ oc 1/ ^Jl^- Therefore, we conclude that this physical process actually cannot happen 
when the cutoff is taken to infinite. 

In summary, this section has discussed the solution for the RPA equations in the context 
of 6 ^ 0+. The excitations are interpreted as quantum particles and the RPA equations are 
treated as a coupled channel scattering problem. The system of equations allows analytical 
solutions and we have obtained the amplitudes of oscillation which are wavefunction of one 
and two-particles. In particular, we have shown that at one-particle channel the effect of 
coupling is to reduce the effective mass of the particles and two-particle channel leads to 
infinitesimal scattering waves for large cutoff. Furthermore, the inelastic processes, cases c. 
and d., cannot occur when the renormalized field theory is considered. 

VI. Discussion 

From the previous sections and I we have seen that the Gaussian variational approxi- 
mation indicates two distinct nontrivial version of (0^)3+1 theories. They are characterized 

22 



by the bare coupling to be infinitesimal of the form b ^^L^^ in the continuum limit. Al- 



though classical intuition and perturbation calculation indicate that the theory with 6 — > 0~ 
is unstable, but preceding analysis and I reveal that it has a well defined symmetric phase. 
In particular the renormalized version of the theory allows a stable vacuum and finite scat- 
tering amplitude [1] [4]. The other one, b 0+, shows a renormalized broken phase solution, 
b(f'^ = 3m^, which is degenerate with the symmetric phase. In section III we have discussed 
the conditions of stability for these vacuums by studying the eigenvalues of the stability ma- 
trix. An alternative approach to investigate this problem is through analyses of the positions 
of the RPA modes. The method is more general because it includes inhomogeneous degree 
of freedom. In the following discussion we will consider the cases of symmetric and broken 
solution separately. 



where we have introduced the covariant variable s = fi^ — [7]. Thus, this channel yields 
a simple free-particle spectrum. The Eq.(5.4), on the other hand, describes the nontrivial 
sector of this vacuum where one has excitations of two particles. 

This case was explored in I and here we will summarize the main results. In order to get 



a. Symmetrical Phase 



a-1. Massive Theory 



In this case 50 decouples from 6p and ( |5.3| ) becomes 
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the spectrum we write the scattering matrix for this process as 

1 



Ti 



p, p', P;fi) = h{p, P) 2 Hp', P) . (6.2) 



--ip{n) 



In appendix B of I we have performed the calculation of I{Q). In terms of covariant variable 
s the result is 



Hs) = -L^ - T^Jis) - 9{s - 4m^)T^V , (6.3) 

lovr^ loTT V s 



where, = l/Svr^ log(2Ap/em) and 



f{s) = 2 + W log^^^- — , s>4m^ (6.4) 



/(.) = 2-2W^tan-M-^ < 5 < 4m^ (6.5) 



„ , , Is — Am? , — 4m2 — ^/s 

f{s) = 2 + J / T-2^ -p 6.6 

V s Vs - 4m2 + y/s 

(cf. fig. 3 of I). Depending on the value of s the system can have different dynamical behavior. 
When s < the system is unstable and for s > 4m^ it has a continuum spectrum. In the 
interval of < s < 4m^ the system may present stable bound states if one finds such as 

A^s,) = l + L^ + ^ = 0. (6.7) 

In this interval, < f{s) < 2, this equation has solution if, only if 6 — ^ In particular 

one can choose the renormalization condition used for the static equations, i.e., ( |2.12D - 



( p.l3| ). In this way, we can show that this theory will result a stable bound state when the 



renormalized coupling constant gm < — I/Stt^ (cf. section V of I). 
a-2. Massless Theory 
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In section II we have shown that with the renormahzation condition a + 6A/167r^ = 
the system has a massless solution. In particular the broken phase is degenerate with the 
symmetric phase. In this case it is still convenient to use fh as the mass scale defined by the 
relation ( p.21| ). Using these ingredients we get 



A+ , < = T - I{s) = log ^ + -— log ^ . 6.8 

In this expression a labels the two theories as we have discussed before. 
Case 1: a = 1 (6 ^ 0") 
The equation ( |6.8|) becomes 

X 1 . 2Ao 1 , A^ 1 , e^m^ ^ ^ 

A+ s < = log ^ + -— log ^ = -— log ^ , 6.9 

oTT^ em 167r^ |s| 167?"^ |s| 

where we have used the relation A^ = (2Ap)^. The system is stable if A"^(s < 0) 7^ for any 
s. However, from (|6.8|) one can always find Si such as A"'"(sj) = 0. Therefore the vacuum is 
unstable when 6 — 0~. 

Case 2: a = -2 (6 ^ 0+) 

In this case the broken phase collapses to the symmetric vacuum. The stability analysis 
is analogous to the previous one: 

X 2 , 2Ap 1 , A, 1 A^ ^ ^ 

^^ 0^ < = ^ log ^ + — ^ log ^ = — ^ log . 6.10 

Stt^ em 167r^ \s\ 167r"^ e*m* s 



The solution for A+(s) = is 



(6.11) 
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This solution is unphysical and has to be ehminated. Therefore, the system is stable in this 
case. 

Thus, the Gaussian variational approximation indicates that (0'')3+i with b ^ 0~ might 
be a nontrivial theory. It has a well defined potential for the sector of (y? = 0, m) of the 
variational space. For certain range of renormalized coupling constant the theory allows 
stable vacuum. Its excitations yield one (free) and two-particle modes. The two-meson 
equation leads to a single bound state and the scattering amplitude in the continuum. We 
also find a massless solution with the renormalization condition a + 6A^/167r^ = 0. The 
RPA analyses indicate that the solution is: i)unstable if 6 — > 0^, which suggest that the true 
vacuum is not homogeneous; ii) stable for theory with 6^0"'". 

b. Broken Phase 

In this case we have the one- and two-particle wavefunction coupled through of bip. In 
order to discuss the RPA equations (|5.10| ) is rewritten in terms of covariant variable s as 

In the interval of < s < 4m^, I^{s) is real and the system may present a stable state if 
one finds a solution > such as 



1^ 3/+(., 



2/6-/+(sJ 
3 



(6.13) 



6(L„ + /(s J/167r2) 

where we have used ( p.3|) . Using now (2.12) with a = —2, we arrive at the foUwing equations 
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for Sg-. 



m 



247r2L, 



1 , ^, 
— + log— + / 



(6.14) 



Thus, one may have a solution for and it is stable if the bracket is positive. From 

fig. 3 of I one can see that f{s) goes to zero for small s. Therefore, we might neglect the last 
term as first approximation to yield 



— 2 

m 



247r2L,, 



+ log 

9fi 



— 2 

m 



(6.15) 



This is precisely the eigenvalue found in Sec. Ill (cf.Eq.(3.9)). 

When s > Affi^ I^{s) has an imaginary part and ( |6.15| ) becomes, after the renormaliza- 



tion. 



— 2 

m 



— 2 

m 



247r2L, 



■(— + log— + /(s) + 



4m2 



Hence, we have continuum spectrum in this case. 



0. 



(6.16) 



It is still illustrative to see the above results from the elastic channel, i.e., /3 — >■ /3, where 
physics is described by the integral equation ( p.l9| ). Note that the potential term is also 



separable and the scattering T matrix can be obtained as usual [8]: 



T(p, p', P;fi) = /i(p, P)S(P;fi)/i(p', P) 



(6.17) 



where S is given by ( |5.22|) . In terms of s it can be rewritten as 



3m 



s — m?' A+fs 



^ s — ' 



(6.18) 



From this viewpoint several dynamical behavior is described by denominator of the scattering 



matrix T, i.e., A+(s) = 0. In particular, the solution found in the equation (|6.15|) can be 
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seen as a bound state of two massless particles. It is stable for < 0. Another intereting 
point is to obsever S(s) when s > Afn^: 



In this range of energy the spectrum is continuum, but the scattering amplitude vanishes 
as In other words, the particles effectively do not interact, but in a very non-trivial 

way. 

In summary, we have obtained RPA equations for 0^ field theory by considering near 
equilibrium dynamics about the critical points of the gaussian parameter space. A simple 
analytical solution can be obtained and it is a nonperturbative method to investigate one 
and two-meson physics. Using this framework we have investigate the stability of the theory 
and shown that in the continuum limit two distinct version of {(p^)2,+i arc viable. One has 
6 — > 0~ and its renormalized theory is stable at the symmetric phase. The other, 6 — > 0+, 
indicates spontaneous symmetric breaking. However, the excitation fields effectively do 
not interact leading to infinitesimal scattering amplitude [9]. Finally, we comment that the 
techniques developed here are general and can be readily extented to other relativistic field 
models [10] and low energy many-body problems such as Boson condensation [11]. 




8 + 2w? 



(6.19) 
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